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We consider two non-linear generalizations of fractal interpolating functions generated from it- 
erated function systems. The first corresponds to fitting data using a 7f th -order polynomial, while 
the second relates to the freedom of adding certain arbitrary functions. An escape-time algorithm 
that can be used for such systems to generate fractal images like those associated with Julia or 
Mandelbrot sets is also described. 



I. INTRODUCTION 



An Iterated Function System (IFS) can be used to construct a fractal interpolating function for a given set of data 
I 0] . The simplest such system defines an IFS 



t ' \ fa n 0\ft\,fe 



c n J \ x J \ /, 



(1) 



with coefficients a n , c n , e„, and /„ determined from discrete data points (U, Xi), i = 0, 1, . . . , N. Such an IFS inter- 
polates the data set in the sense that, under certain assumptions on the coefficients [Q, the attractor of the IFS is a 
graph that passes through the data points. In this particular case, the IFS can be written as 



t , = (t - t ) t + (t- t N ) 



(tN — to) (to — tN) 

(f ~ t n -l) - t n ) 



(2) 



which shows that a linear (in t) interpolating function between the points (t n _i, x n -i) and (t n ,x n ) is used. 

Various generalizations of fractal interpolating functions have been given, including those for higher dimensional 
functions, the use of hidden variables, and extensions to certain non-linear distortions |3|, |], ^|, ^|, In this note 
we describe a generalization whereby the transformation incorporates a JC th -order polynomial interpolation between 
adjacent points. We also discuss certain classes of non-linear functions that can arise in such interpolating functions, 
and show how such functions can, with the use of a particular escape-time algorithm, be used to generate certain 
fractal images. 

The paper is organized as follows. In Section || we describe simple linear fractal interpolating functions, and discuss 



how particular non-linear functions can arise. Section [II generalizes these considerations to if —order interpolating 



functions. Section IV describes a certain escape-time algorithm which may be used for these systems to generate 



fractal images like those associated with Mandelbrot or Julia sets. Section M contains some brief conclusions 



II. LINEAR INTERPOLATING FUNCTIONS 



We first describe how a standard linear fractal interpolating function is constructed. Suppose we have data points 
(ti, Xi), i — . . . N, describing a function x(t). Consider the IFS 



Imposing the conditions, for n = 1, 2, . . . , JV, 



to \ _ ( tn-l 

Xo I \ Xn-l 
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leads to determination of the coefficients as 



tn tn—1 

O-n = 

tN — to 

tjM — t 
Xn %n— 1 
tN — to 

X n -itN — X n t . . 

In = " T (5) 

tjv — to 



The transformation can then be written as 



W (t\ - t' - ^ ~ to ) t 4- ( t-tjv ) f 
W n {t)=t — — — t n + — — t„_i 

{tN — to) [Zq — tjvj 

W n (x) = X = j X n + -H —r X n -l (6) 

[tn tn— x) [tn—1 tn) 

Thus, W n (x) = x 1 is determined by a linear (in t) interpolating function constructed between the points (t n _i,a; n _i) 
and (t n ,x n ). 

A generalization of this type of fractal interpolating function can be found by considering an IFS of the form 

W n (t) = a n t + e n 

W n (x) = c n t + f n + g n (x) (7) 

where g n (x) is, at this stage, an arbitrary function. Imposing the conditions (Q) leads to determination of the 
coefficients as 

tn t n —\ 



tN — to 
tn-ltN — t n to 
tN — to 

^ _ X n - Xn-l _ 9n(x N ) - g n (xp) 

ijv — to tN — t 

j. _ Xn-ltN — X n to 9n{xo)tN ~ 9n{xN)to /gs 

tN — to tN — t Q 



The transformation can then be written as 



where 



and the identity 



W n (x) ee x' = £ x n + — ^j— x n -i (9) 



*n = + g n (x) - g n {x N ) 

Xn-l = Xn-l + 9n(x) - g n (xo) (10) 



f-^U/-^ =1 (11) 



for arbitrary has been used. 



3 



III. QUADRATIC AND HIGHER ORDER INTERPOLATING FUNCTIONS 

The interpolating function of the last section used a linear (in i) approximation between adjacent points. In this 
section we indicate how a quadratic approximation may be constructed; the generalization to an arbitrary if^-order 
polynomial approximation will be straightforward. Let us consider a transformation of the form 

W n (t) = a n t + e n 

W n (x) = c n t + d n t 2 + f n (12) 



and impose the conditions, for n = 2, 3, . . . ,N, 

Wn 



to \ _ / t n -2 
XO J \ X n -2 

tm \ / tn—1 

Xrn J V X>n — 1 



XN I V x 



Wn ( I* ) = ( t n ) (13) 



The point t m is determined as 



(tn-l-tn-2) , ftn-l _ *n) 

tm = "77 7 r~ tN + 77 7~T t{) \ lA > 

V-n — E«-2j V"n-2 — In ) 



with corresponding point x m . The coefficients of the IFS are determined as 

tn — t n -2 

a n = 

tff — to 

_ t^n-2 — t()t n 
tN — to 

_ x n(tp — t m ) + X n -i(t N — tp) + X n ^2(t m — tjy) 
(tN — to)(tN — t m )(t m — to) 
j _ X n (t m — tp) + X n -i(tp — t N ) + Z n _ 2 (tjy — t m ) 

(tN — to)(tN — t m )(t m — t ) 
j. _ X n t m to(tm — to) + X n -\tNto(to ~ tN) + X n -2tNt m (tN — t m ) 
(tN - t )(tN - t m )(t m - tp) 

With this, the transformation can be written as 

w (a ft - to) . , ft ~ *n) , 

W n {t)=t - — — t n + — — — t n -2 

(tN — to) (t — tN) 

tit ( \ — I _ ft' ~ t n -l)(t' — t n -2) ft' — t n )(t' — tn-2) 

(t n — t n -i){t n — t rl -2) (t n -i — t n )(t n -i — t. n -2) 

(f -tn-l)(t' -t n ) 

X n -2 (16) 



(tn-2 — t n -l)(t n -2 — t n ) 



which thus uses a quadratic (in t') interpolating function between the points (t n ,x n ), (t n -\, x n -i), and (t n -2,x n -2)- 
As in the previous section, including an arbitrary function g n (x) in the IFS transformation via 

W n (t) = a n t + e n 

W n (x) = c n t + dj 2 + f n + g n (x) (17) 

is straightforward. The conditions (|l3| ) leads to determination of the point t m of Eq. (|l4| ) as before, together with 
the accompanying point x rn . The transformation itself can be written as 

W (t\-t> {t - to) t 1 ( *~ tjv) t 

W n (t)=t - — TTt n + — — -rr t n -2 

(tN — t ) (to — tN) 

TJT { \ _ I _ (t' — t n -l)(t' — tn-2) _ (t — t n ) (t — tn-2) 

(t n — t n -i)(t n — tn-2) (t n -i — t n )(t n -i — tn-2) 

X n -2 (18) 



(tn-2 — t n -l)(t n -2 — t n ) 
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where 



x„ = x„ +g n (x) - g n (x N ) 

X n — \ — X n —i ~t~ Q n {x) ~ Qni&rn) 

x„-2 = x n - 2 + g n (x) - g n (x ) (19) 



and the identity 

(* / -tn-l)(* / -*n-2) (*' -t n )(f-t n - 2 ) {f -t n -l)(t' -t n ) 



(20) 



(t n — t n -l)(t n — t n - 2 ) (tn-l — t n )(t n -l — t n ^2) (fn-2 — t n -l)(t n -2 ~ Ui) 

for arbitrary t' has been used. 

From these considerations, the pattern to constructing a iif th -order fractal interpolating function is apparent. Start 
with a transformation of the form 

W„(t) = a n t + e n 

Wn {x) = + B$t + B^t 2 + ... + Bi K h K (21) 
and impose the conditions, for n = K, K + 1, . . . , N, 

w n 
w n 



to \ _ I t n ^K 

%0 J ~ \ X n - K 

tml \ _ I t n —K+\ 

Xml ) V Xn-K+1 



t m 2 \ ( t n -K+2 

X m 2 I \ X n -K+2 



t]Sl \ I t r . 



»'■(*: ) = (z) < 22 > 

The K — 1 intermediate points t m j, with j = 1,2,. . . , K — 1, are determined as 

_ (tn-K+j - tn-x) (tn-K+j - tn) . . 

t m j - 7- " \ t N + — — — t {16) 

(t n — t n -K ) (tn-K — t n ) 

along with the corresponding x m j points. The resulting transformation will be of the form given by Lagrange's formula 
for a if th -order polynomial interpolating function constructed from K + 1 points: 

W n \t)=t - -r- rT% + 7T —:t n -K 

(in — to) (to — tN) 

w r s _ I _ jt' ~ t n -l){t' - t n _ 2 ) ■ ■ ■ [f ~ tn-K) 
W >fi ( X i — X — , , , x, v Xfi I 

(tn — t n —\)(t n — tn.— 2) ' ' ' (t n — t n -K) 

(t — t n )(t — t n _ 2 ) 4 4 ■ (t — t n _ K ) 

r X„_i + . . . + 

-K) 

.A 

x„- K (24) 



(t n -l ~ tn)(t n -l — t n -2) ■ ■ ■ {tn-1 — tn-K) 
[t — t n )(t — t n -l) ■ • ■ (t — tn-K+l) 



(t n -K — t n )(t n -K — t n -l) • ' • {tn-K — t n -K+l) 



The inclusion of an arbitrary function Qn(x) in the transformation W n (x) of Eq. (p2|), as was done for the linear 
and quadratic transformations of Eqs. (0) and ([l?]) respectively, is straightforward. As might be expected, the use 
of these higher-order interpolating functions can increase the accuracy of the interpolation significantly, at least for 
smooth functions - some informal tests on known functions suggest an improvement of almost an order of magnitude 
in general in using a quadratic interpolating function over a linear one. Of course, as for polynomial interpolation, 
there is a limit to the net gain in employing a higher-order interpolating function. 
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IV. ESCAPE-TIME ALGORITHM 



Assuming that the corresponding IFS transformation is contractive, so that the distance d(t, x) between any two 
points in the range of interest satisfies 

d(W n {t),W n {x)) <s n d(t,x), (25) 



where < s n < 1 is the contractivity factor, graphs of the functions represented by fractal interpolating functions 
can be made by applying the standard random iteration algorithm to the IFS: 

• initialize (t, x) to a point in the interval of interest 

• for a set number of iterations 

- randomly select a transformation W n (t, x) 

- plot (t',x') = W n (t,x) 

- set (t, as) = (t',x') 

• end for 



Alternatively, one can relate an IFS W n (t, x) to a shift dynamical system f(t,x), and on this system perform an 
escape time algorithm to generate an image fy] . In this section we describe an algorithm for generating fractal images 
like those for Julia or Mandelbrot sets from IFS interpolating functions. 

Suppose we have an IFS transformation W n (t,x), generated by some data points {Xi,yi), £ = 0,1,. which 
includes a non-linear function g n (x), as was done for the linear and quadratic transformations of Eqs. (Q) and ( |l7j ) 
respectively. We now continue the real variable x of this transformation to complex values: x — > z = (zr, zj), so that 
the transformation W n (t,z) is defined on the complex plane. We can then, in analogy with the algorithm used for 
Julia sets, define the following escape-time algorithm to generate a fractal pattern: 

• for each pixel in a region of interest 

— initialize t 

— initialize z — {zr,Zi) to the pixel coordinates 

— for n = 0, l,...,N 

* calculate (t',z r ) = W n (t,z) 

* break if \J z'^ + z' 2 exceeds a maximum 

* set (i, z) = (<', z') 

— end for 

— plot the pixel 

• end for 



where the pixel is plotted using a coloring algorithm based upon, amongst perhaps other factors, the number of 
iterations attained when the break condition was met S. 

The preceding can be interpreted as follows. A general A^ th -order fractal interpolating IFS 

W n (t) = a n t + e n 

W n (x) = 4°) + B^t + B^t 2 + ... + Bi K h K + g n (x), (26) 

with the coefficients a n , e n , B^ , B^ , ■ • ■ , B^ determined from the data (t », Xi), i = 0, 1, . . . , N, can be viewed, with 
the continuation defining a complex map 

z n+1 = B^ + BgK + B^tl + ... + Bi K h« + g n (z n ) (27) 



for n = 0, 1, . . . , N. The escape-time algorithm described above is then just the standard one used for Julia sets of 
complex maps. 
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The arbitrariness of the function g n (x) and the data set (tj, Xj) used to fix the IFS leads to a wide variety of possible 
fractal images generated in this way. An interesting class of functions g n (x) to consider in this context are those for 
which, when continued to the complex plane x — > z = (zr, Zi), lead to a map having a fixed point z\ = 0: 

Zj = = lmg n (z R , z*j) (28) 

In such a case one could augment the usual condition of the escape-time algorithm to cease iteration: \/ z 2 R + z| > A, 
where A is some suitably large value, to also cease iteration when \zj\ < e, where e is a suitably small value. The 
coloring algorithm used to plot a pixel, which depends on the number of iterations attained when this break-out 
condition was met (if at all), will then lead to structure in the region where the break-out condition on the magnitude 
of z is not met. 

We give two examples of fractal images generated this way for the choice g n {x) = —0.4a; + 0.5x 2 + 0.2a: 3 , with the 
data generated from the logistic map x n +i = 3Ax n (l — x n ), with n = 0, 1, . . . 60. The first one, appearing in Fig. fiL 
corresponds to the generalization (|7|) of a linear (in t) fractal interpolating function, while the second image of Fig. g 
corresponds to the generalization (TPa) of a quadratic (in t) interpolating function. A coloring algorithm that simply 
mapped a color to the number of iterations attained when the break-out condition became satisfied was used in both 
cases. 




FIG. 1: Fractal image from a f-linear interpolating function 
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FIG. 2: Fractal image from a t-quadratic interpolating function 

These figures illustrate, in the interior of the fractal object, the richer structure arising from the quadratic over 
the linear interpolation function. In this region the break-out condition \zt\ < e is satisfied, which numerically for 
e ~ 10~ 5 is attained after a relatively small number (10-30) of iterations. 

V. CONCLUSIONS 

We have considered two non-linear generalizations of fractal interpolating functions constructed from iterated 
function systems. One - using a i4T th -order interpolating polynomial - can potentially improve the accuracy of fractal 
interpolating functions. The other generalization - the use of certain arbitrary functions in the IFS - can, together 
with an appropriate escape-time algorithm, generate fractal images. This last point is of interest as, first of all, there 
is a rich variety of such images possible due to the arbitrariness of the functions used, and secondly, it shows how 
fractal images as normally associated with Julia or Mandelbrot sets can also be associated with discrete data sets. 
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